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Abstract 

The motion of a particle in a spatially harmonic magnetic field is a basic problem involved, for 
example, in the mechanism of formation of a collisionless shock. In such settings, it is generally 
reasoned that particles entering a Weibel generated turbulence are trapped inside it, provided their 
Larmor radius in the peak field is smaller than the field coherence length. The goal of this work 
is to put this heuristic conclusion on firm ground by studying, both analytically and numerically, 
such motion. A toy model is analyzed, consisting of a relativistic particle entering a region of 
space occupied by a spatially harmonic field. The particle penetrates the magnetic structure in a 
direction aligned with the magnetic filaments. Although the conclusions are not trivial, the main 
result is confirmed. 
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FIG. 1: Setup considered. The field is 0 for z <0 and Bq sin(kx)e y for z > 0. Particles arrive in 
the field region with x(f = 0) = (xo, 0,0) and x(i = 0) = (0,0, vo)- 

INTRODUCTION 


The filamentation instability, triggered when two counter-streaming plasmas collide, 
grows a nearly spatially harmonic magnetic field [1, 2], When it reaches saturation, it 
leaves an array of magnetized filaments which start merging ^-7]. If plasma keeps flowing 
in the region occupied by this field, the question of the particles’ trajectories in it comes 
into play. This issue is specially relevant to the mechanism of collisionless shocks formation. 
In this respect, Particle-In-Cell simulations have shown that the shock starts forming when 
the peak magnetic field, together with its coherence length (1 /k), are large enough to stop 
the incoming flow in the region where the instability grew 8j, 9]. An intuitive case, which 
proves correct, has been made that for this the occur, the field coherence length, i.e. the size 
of the magnetic filaments, must be larger than the Larrnor radius of the incoming particle 

in the peak field HQ. 

The goal of this work is to put this reasoning on firm ground, by providing a rigourous 
description of the particles’ trajectories in a spatially harmonic field. To do so, we model 
the process as pictured on Figure [0 Note that the very nature of the problem forbids us to 
use any slowly varying field concept, like the adiabatic invariant for example 12]. 

We consider a particle of mass m and charge q travelling through an inhomogeneous 
magnetic field. The field B reads 0 for z < 0 and Bo sin(kx)e y for z > 0. At t — 0, 
the position of the particle is x(i = 0) = (x 0 ,0,0) and its velocity x(i = 0) = (0,0,u 0 ), 
with v 0 > 0. Note that we conduct a test-particle type approach of the same kind that 
the one conducted in the heuristic reasoning described earlier. Collective effects could be 
investigated in further works. Further caveats of the model are listed in the conclusion. 
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Since the particle is subjected to a magnetic field only, its Lorentz factor 7 is constant. 
The equation of motion for z > 0 then reads, 

x 

myx = q— x B. (1) 

c 

Setting, 

ojb = ———, x— >X./k, t^r/cuB- (2) 

'ymc 

Equation (JT() reads for Z > 0, 

X = -Z sinX 
Y = 0 

Z = X sin X (3) 

with initial conditions, 

X(t = 0) = (kx 0 , 0, 0), 

X(r = 0) = (o,0,^£). (4) 

The initial dimensionless velocity Z 0 = kv o/ojb > 0, is simply the ratio of the Larrnor radius 
to the filament size. 

The second equation, Y = 0, coupled with the initial conditions above, simply means 
that the motion lies in the (. X , Z) plane. The problem is bi-dimensional. The differential 
system under scrutiny is eventually, 

X = -ZsmXH(Z), (5) 

Z = Xsin XH(Z), ( 6 ) 

where H(Z) is the step function, H(Z) = 0 for Z < 0 and H(Z ) = 1 for Z > 0. The initial 
conditions are defined from (@| as, 

X Q = kx 0 , 

Z 0 = kv 0 /uJ B ■ (7) 

Before we move to some analytical results, let us plot a few trajectories by numerically 
solving our system. To this extent, Eqs. (1515171) have been solved using the NDSolve 
function of the Mathematica software. Figure [7] shows how, depending on their initial 
conditions, some trajectories return to Z < 0, while others proceed to Z — + 00 . Our goal 
is to analyze the ultimate evolution of the trajectories in terms of their initial conditions. 
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FIG. 2: (Color online) Numerical solutions of the system (13141) . The red and the purple trajectories 
both have Xq = 2. But the red one has Z$ = 0.1 while the purple one has Zq = 0.9. The beige 
trajectory starts with the same velocity than the purple one, Zq = 0.9, but from Xq = 1. The red 
trajectory returns to Z < 0 while the others proceed to Z = +oo. Starting from Xq = n ± 0.1 
respectively, the green and blue trajectories are bounded in the X direction by the inequality Q. 
The horizontal line at X = it is never crossed (see text). 



FIG. 3: (Color online) Trajectories are bounded by Eq. d9j) within intervals of width < 27r. 

ANALYTICAL RESULTS 

X-range and -periodicity 

We here establish that the motion is always bounded in the X direction, which will 
eventually allow to restrict our study to X 0 e [0,7r]. Integrating (JH]) from 0 to r gives 

Z (r) — Z 0 = cos X 0 — cos X. (8) 
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Since the velocity is constant in the magnetic held and Ao = 0, we know that Z(t) < Zq. 
Therefore, and as a consequence of the equation above, cos A 0 — cos A < 0, i.e, 

cos A" > cosA 0 . (9) 

Depending on the sign of sin A 0 , the particle is initially dehected upward (if sin A 0 < 0), or 
downward (if sin Ao > 0). Figure [3] illustrates how these constraints result in a confinement 
of the trajectory within intervals of width < 2n. As a consequence, we can restrict our study 
to Ao G [—7T, 7r]. In addition, the system (I5ll6l) is clearly invariant when changing X —> —A, 
so that it is enough to investigate the interval A 0 G [0,7r]. 

Sufficient condition for linv^oo Z(t) = oo 

Integrating now ([8]) from 0 to r and accounting for Zq = 0, we find, 

Z(t) = f (Z 0 + cos A q — cos X)dr. (10) 

Jo 

From — 1 < — cos A < 1, we deduce the following inequality for the integrand, 

Zq -f- cos A o — 1 A Zq T cos A o — cos A A Zq T COS Aq T 1. (11) 


Therefore, 

• If Zq + cos A 0 — 1 > a for any a > 0, then Z(t ) > ar and lim r _ > . 00 Z(r ) = oo. The 
particle streams through the magnetized region. 

• If Zq + cosA 0 + 1 < —a for any a > 0, then Z(r) < —ar and lim^oo Z(r) = —oo. 
The particle bounces back on the magnetized region. 

Note that in the second case, our analysis is only valid as long as Z > 0. Yet, as soon as 
the particle with an initial Zq > 0 crosses back the frontier Z = 0 from Z > 0, it finds itself 
in a field-free region where it won’t be deflected again toward a positive value of Z. Hence, 
our conclusion remains unchanged: the particle returns back to Z = —oo. 

The first conclusion is already a partial confirmation of the intuitive reasoning that the 
threshold for bouncing back is Z 0 ~ 1. Averaging over the possible values of A" 0 , we find 
particles stream through the magnetized region if Zq > 1, i.e. Vq/lub > k -1 . This is precisely 
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equivalent to stating that the filaments are too small compared to the Larmor radius of the 
incoming particle in the peak held. 

The second conclusion is indeed trivial since Z 0 + cos A 0 + 1 < 0 eventually amounts to 
Z 0 < — cosA 0 — 1 < 0: if the particle is launched from Z 0 = 0 with Z 0 < 0, it necessarily 
goes to Z = —oo, or stay at Z = 0 (in the case Z 0 — 0). 


Linearization for sinX <C 1 


A more accurate treatment of the system (151161) can be performed if sinX 0 ~ 0. Let us 
start assuming that sinXo ~ 0 implies sinA(r) ~ 0, Vr > 0, and check the assumption 
afterward. If X 0 ~ 0, we’ll then write sin X ~ X. If Ao = n, we can replace X —> X — n and 
analyze the system considering sin A ~ —A. The system of linearized equations eventually 
reads, 

A = -eZX, (12) 

Z = eXX, (13) 


with e = 1 if A 0 ~ 0, e = — 1 if A 0 ~ i r, and initial conditions (|4|) . Integrating Eq. (TT3l) 
gives, 

Z=Z 0 + e -X 2 - e -Xl (14) 

Note that this equation is nothing but the linearization of Eq. (JH]) for sin A" <C 1. Inserting 
this expression of Z in Eq. (TT2|) and neglecting second order A terms and superiors (including 
Aq), we find, 

A + oj 2 X = 0, u> 2 = eZ 0 => X(t) = A 0 cosh(io;r). (15) 


For A"o ~ 7r, e = — 1 and iui is real. Here, trajectories quickly depart from the region 
where sin A" ~ —A", and linearization breaks down. The trajectory (A = 7 t,Z — Z 0 t ) is 
unstable. But for A 0 ~ 0, e = 1 and A (r) = A 0 cos (car): trajectories remain conhned 
within X G [—A 0 ,A 0 ], and the linear treatment remains valid. Using then Z = A A and 
integrating twice gives, 


Z{r) = (Zo - ^ 


sin(2 cur) 2 
T + —x-A 0 . 


(16) 


4 ) 8u 

The particle is thus repelled from the magnetized region if Z 0 < Ag/4, a slightly more 
stringent condition than Z 0 < 1 — cos A 0 Aq/ 2, near A 0 = 0. 
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FIG. 4: (Color online) Whenever the initial conditions (Zq, Xq) pertain to the shaded area, the 
particle bounces back. 

NUMERICAL RESULTS 

In order to derive, for any A 0 , the sufficient and necessary condition for particles to bounce 
back against the magnetized region, we numerically explore the phase space (Zq,Xq). For 
any given A 0 € [0,7r], we scan the values of Zq yielding the particle to bounce. 

The result of this numerical exploration is summarized in Figure |4j whenever the initial 
conditions (Z 0 , X 0 ) pertain to the shaded area, the particle bounces back. The full picture 
can be completed by symmetry around X 0 = 7 r, and X 0 — 0. 

For Xq ~ 0, we numerically find the bouncing threshold at Zq ~ 0.34A"q . The discrepancy 
with the Zq = Aq/4 condition derived in the linear regime is explained as follow: some 
forward trajectories, quite like the purple one on Fig. [21 come close to exit the Z > 0 
region by their beginning. Even if they would go to Z — +oo in a space hlled-up with the 
harmonic held, they leave the magnetized region early in time, simply because they wander 
into the field-free zone Z < 0 during their first oscillation. Trajectories with A 0 ~ 0 are 
forward (backward) for Zq > (<)A"q/ 4. But they wander into the held free zone as soon as 
Zq < 0.34Ag. 

This factor 0.34 can be approximated analytically in the following way. One starts lin¬ 
earizing Eq. (TT6j) up to the 3rd order in the time parameter r. Equaling the result to 0 gives 
the moment r 0 when the particle, starting from Z = 0, will reach Z = 0 again. The result 
is r 0 = V^Zq/uXq Q- Then, if Z(r 0 ) < 0, the particle will dehnitively escape toward the 
Z < 0 region. The time derivative of Eq. (1T6T) now readily gives Z(r). Setting r = tq in the 
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FIG. 5: Portion of a spatially uniform incoming beam bounced back against the magnetized region, 
in terms of Zq. 

result and developing it up to the 3rd order in Z a gives Z(r 0 ) < 0 O Z 0 < |Xq, which is 
reasonably close to Z 0 < 0.34Xq. 

For X 0 ~ 7 r the shaded area displays a “hole” between Z 0 = 0 and 1. We here deal with 
trajectories like the blue or the green ones on Fig. [2] As stated earlier, the linear analysis 
fails here because trajectories quickly departs from A" ~ it. But the numerical investigation 
shows that these trajectories starting from A 0 ~ 7r have a common point with those starting 
near Xo: because the held is zero a both locations, particles starting exactly from Xo = n 
or 0 definitely stream through. When starting slightly aside from these values, it only takes 
them a small initial velocity to override the held and move forward. 

This gap in the bouncing region forbids to average over Xo and derive an averaged 
bouncing, or streaming through, condition. We can instead assume an homogenous incoming 
beam covering all values of X 0 , and numerically derive from Fig. [T] the portion of it which 
is bounced back, in terms of Zq. The result is pictured on figure [5] For zero initial velocity, 
100% is bounced back. The percentage then decreases until it reaches 0 for Z 0 ~ 1.09. We 
could not find a way to approximate this number analytically. Macroscopically, we thus 
recover a condition for the beam to stream through the magnetized region, of the form 
Z G > 1, that is, k _1 < vq/cob- 



















CONCLUSION 


When dealing with more realistic settings like the ones arising from Particle-In-Cells 
simulations, the setup is far from being so simple. In this respect, the model should be 
refined in at least 3 ways. 


1. The hlamentation instability does not grow a single mode, but a continuum of modes 
instead. Although the fastest growing one in the linear phase definitely governs the 
end picture, some smaller and larger wavelengths also grow at a slightly slower rate 
Q, so that the resulting held is more involved than the present single harmonic one. 


2. The transition from the held-free region to the magnetized one is not clear cut, but 
displays an involved gradient instead Q. ,n addition, the existence of electric fields 


generated by charge separation has been demonstrated in this transition layer 15]. 
Finally, still in the transition region, the longitudinal wavelength is comparable to the 
transverse wavelength so that the semi-infinite filaments model fails. 


3. It is assumed that B = 0 for 2 < 0. I 11 a realistic shock configuration, this might be 
identified with the upstream region. However, the upstream region is not unmagne¬ 
tized, but it has magnetic filaments that are roughly stationary in the fluid frame. This 
implies that, for a particle moving with the E x B drift, the force will vanish, similarly 
to the case of a region with zero electric and magnetic fields (as assumed in this work). 
However, the situation will be different for particles reflected back upstream by the 
shock. In the present model, they will continue towards upstream infinity, whereas in 
a more realistic shock they would be advected back towards the shock. 


The present model seems therefore more appropriate to the early stages of formation of 
a collisionless unmagnetized shock. Here, particles with Z 0 > 1, i.e. Vq/lub > /c -1 , mostly 
stream through the magnetized region 8]. But those with vq/ujb < k are allowed to pass 
through the interface, before they are trapped in the magnetized region (instead of bouncing 
back). This is the onset of the collisionless shock formation. At any rate, we have checked 
that Vq/ojb ~ k~ v is the threshold for two completely different qualitative behaviors. The 
incoming flow gets trapped in the Weibel generated turbulence, if the Larmor radius of the 
particles in the peak held is smaller than the size of the filaments. 
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